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We study analytically a simple random walk model on a one-dimensional lattice, where at each time step the 
walker resets to the maximum of the already visited positions (to the rightmost visited site) with a probability r, 
and with probability (1 — r), it undergoes symmetric random walk, i.e., it hops to one of its neighboring sites, 
with equal probability (1 — r)/2. For r = 0, it reduces to a standard random walk whose typical distance grows 
as i/n for large n. In presence of a nonzero resetting rate 0 < r< 1, we find that both the average maximum and 
the average position grow ballistically for large n, with a common speed v(r). Moreover, the fluctuations around 
their respective averages grow diffusively, again with the same diffusion coefficient D(r). We compute v(r) and 
D(r) explicitly. We also show that the probability distribution of the difference between the maximum and the 
location of the walker, becomes stationary as n —>■ oo. However, the approach to this stationary distribution is 
accompanied by a dynamical phase transition, characterized by a weakly singular large deviation function. We 
also show that r = 0 is a special ‘critical’ point, for which the growth laws are different from the r —> 0 case 
and we calculate the exact crossover functions that interpolate between the critical (r = 0) and the off-critical 
{r —> 0) behavior for finite but large n. 

PACS numbers: 05.40.-a, 02.50.-r,87.23.Ge 


I. INTRODUCTION 

Search problems appear in diverse contexts and there 
have been a recent surge of interests in the physics community 
in these problems 0. Search strategies may be either system¬ 
atic or random. In the systematic strategies, the searcher uses 
deterministic rules (e.g., ‘lawnmower’) to find a target. On the 
other hand, the random search mechanism typically involves 
two kinds of moves: local steps when the searcher looks for a 
target, and long-range moves during which the searcher does 
not look for the target but relocates itself to a different terri¬ 
tory. The slow search phase is typically modelled by a diffu¬ 
sion or a random walk. The long-range moves may be mod¬ 
elled depending on the specific application IH. 

A particularly simple long-range strategy consists of ‘reset¬ 
ting’ the searcher to a fixed location (say to the initial starting 
point) with a finite probability/rate. The rational behind this 
strategy is that if one does not succeed in finding the target 
via short-range diffusion, it is better to ‘restart’ the process, 
rather than continuing on the short-range moves. The effect 
of such stochastic resetting was first studied by Manrubia and 
Zanette m in the context of multiplicative processes and a 
slightly different version was studied later by Gelenbe Q in 
the context of network theory. Such ‘restart’ strategy also 
plays an important role in randomized search algorithms for 
combinatorial optimization problems ||8] |9l . 

Recently, a very simple model of a Brownian searcher in 
presence of stochastic resetting to its initial position with rate 
r was introduced by Evans and Majumdar ifTOll . In presence 
of a nonzero r, it was shown that at long times, the probabil¬ 
ity distribution of the position of the walker reaches a non¬ 
equilibrium steady state M- The temporal relaxation to this 
steady state was also studied recently ini and an interesting 
dynamical phase transition was found; as time progresses, an 
inner core region around the resetting point reaches the steady 
state, while the region outside the core is still transient. The 


boundaries of the core region grow linearly with time at late 
times lim . 

In presence of resetting with rate r, the mean first-passage 
time to find a target located at the origin, by a searcher starting 
and resetting to xq, was computed exactly m and was found 
to have a minimum at an optimal resetting rate r*, thus mak¬ 
ing the search process efficient in presence of resetting. This 
conclusion holds in all dimensions ifT^ . Also, it was proved 
that this non-equilibrium reset dynamics is more efficient in 
target search compared to an equilibrium Langevin dynamics 
in presence of an external potential leading to the same steady 
state IIT 3 I . 

This simple model of diffusion with stochastic resetting, in 
the single searcher setting, has been generalized in various 
ways. For example, when the target as well as the resetting 
positions are not fixed but drawn from specified probability 
distributions M, in presence of partial detection (or absorp¬ 
tion) of the target by a searcher na, when the searcher per¬ 
forms a continuous-time random walk m or a Levy flight 
instead of a Brownian motion/random walk imiisi, when 
the searcher moves in a bounded domain im or in the pres¬ 
ence of a confining potential ||20| . when the resetting occurs to 
any of the previously visited sites with a rate proportional to 
the number of visits to the site ED etc. Recently the model of 
random walks with resetting has also been used to understand 
the behavior in models of enzymatic reactions in biology fT2\ . 

Going beyond the one particle setting, the effect of the 
resetting mechanism in searching an immobile target has 
also been studied in presence of multiple, but non-interacting 
searchers Eol. More recently, the resetting has been studied 
in spatially extended many-body interacting systems, such as 
for fluctuating interfaces iniiii as well as a class of reaction- 
diffusion models ll24l — in both cases, the natural dynamics of 
the system is stochastically interrupted by resetting it to the 
initial configuration at a nonzero rate r. A nonzero r leads to 
new non-equilibrium stationary states in such extended sys- 
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terns mniia. 


II. THE MODEL AND THE SUMMARY OF MAIN 
RESULTS 


In this paper, we consider a model where the searcher re¬ 
members the maximum location visited so far and the long- 
range move consists in resetting to this current maximum. 
This strategy may be thought as the mixture of the system¬ 
atic search and the random search. In the systematic search 
strategy since each location is visited only once, if a target in 
an already visited place is missed (by an imperfect searcher), 
it is never going to be detected. In the new strategy discussed 
in this paper, the searcher revisits already searched locations 
(with certain probability), but also feels a dynamical bias to¬ 
wards exploring new locations (by resetting to the maximum). 


This model may actually be useful in the context of ani¬ 
mals searching for food. During the foraging period, it is well 
known that an animal typically performs a random walk in 
search of food Il25ll2^ . It is however quite natural for an in¬ 
telligent animal (with memory) to remember the already vis¬ 
ited (explored) sites and thus to have a natural tendency to 
relocate once in a while to the frontier between already ex¬ 
plored and yet unexplored territories, where the probability of 
hnding food may be higher due to the proximity of the un¬ 
explored territory. In a one dimensional setting, the animal, 
besides short-range diffusion, may relocate with a nonzero 
probability to the current maximum or to the current mini¬ 
mum, which together constitute the frontier between explored 
and unexplored territories. In this paper, we consider an even 
more simplified directed version that has the advantage of be¬ 
ing exactly solvable. In our model the animal starting at the 
origin, besides performing short-range standard random walk, 
relocates stochastically with a nonzero probability to only the 
positive side of the frontier, i.e., the farthest visited site so far 
to the right of the origin (i.e., to the maximum). We will see 
that despite the fact that the position of the walker evolves via 
a non-Markovian dynamics (as it remembers the maximum 
position so far in order to relocate), the model allows for an 
exact solution and thus provides interesting insights into this 
long-range search strategy. 


The rest of the paper is organized as follows. In Section [n| 
we introduce the model precisely and summarize our main ex¬ 
act results. In Section III we derive the generating functions 
for the probabilities of the position and the maximum after n 
time steps. We then examine in detail the asymptotic large 
n statistics in the two opposite limits respectively in the next 
two sections: (i) the case without resetting but only with diffu¬ 
sion in Sectionj^and (ii) the case where only resetting occurs 
without any diffusion in Section]^ In Section VI we analyze 
the situation for arbitrary resetting rate. In Section |VII[ we 
analyze the statistics of the maximum and the position in the 
scaling limit r^0,n^°° while keeping the product rn fixed. 
The crossover scaling functions are computed exactly in this 
Section and compared to numerical simulation results. Finally 
we conclude with a summary and some open questions in Sec¬ 
tion Ivm] 


We consider a walker moving on a one-dimensional lattice, 
initially starting from the origin. Here each lattice site should 
be thought of as a ‘region’ which is much larger than the ‘size 
of the searcher’, but much smaller than the whole region. The 
searcher spends some characteristic time T in each region (lat¬ 
tice site), and we consider time steps in units of T and take it 
to be discrete. Let x{n) denote the position of the walker at 
step n and m{n) denote the current position of the maximum 
at step n, i.e. 

m{n) = max [x(0) = 0, x(l), x(2 ),... ,x{n)]. (1) 

The position x{n) evolves with time via the following 
stochastic dynamics. At any given time step n, if the position 
x{n) of the walker is less than the maximum position m{n) 
reached up to that time (i.e, x{n) < m{n) strictly), then in 
the next time step, the position is reset to the maximum 
position with probability r. With the remaining probability 
(1 — r), the walker moves either to the right or to the left 
lattice site, with equal probability (1 — r)/2. On the other 
hand, if x{n) = m{n), then in the next time step, the walker 
moves either to the right or to the left lattice site with equal 
probability 1 /2. The dynamics is precisely dehned by the 
following evolution rules [see Fig.[T](a)]: 


if x(n) < m{n) 


{ (x -I- 1, m) 
(x — 1, m) 
{m,m) 


with probability (1 — r)/2, 
with probability (1 — r)/2, 
with probability r, 


( 2 ) 


and if x(n) = m{n). 


(x,m) — 


(jn — I, m) 


with probability 1/2, 
with probability 1 /2. 


(3) 


Evidently, the evolution of x(n) is non-Markovian by itself, 
since the walker has to remember the maximum position 
reached so far in order to reset. However, the dynamics of 
the pair of stochastic variables {x(n),m(n)} is Markovian in 
the two dimensional (x, m) plane, and this is the key point be¬ 
hind the solvability of the model. Figure [^depicts the motion 
in the (x, m) plane. 

It is useful to summarize our main results. Our main objec¬ 
tive is to compute the statistics of the two random variables 
x{n) and m{n). Let us first recall that in the absence of reset¬ 
ting (r = 0), the walker performs a standard one dimensional 
random walk, for which x(n), converges to a Gaussian random 
variable with zero mean and variance n, for large n. Hence, the 
probability distribution of the position converges, for large n, 

to Px{x,n) —>■ [~-*^/2«] ■ Similarly, for r = 0, the 

distribution of the maximum m(n) > 0, for large n, converges 
to a half-Gaussian: P^[m,n) —>■ exp [—m^/2«] with sup¬ 

port only over m > 0. Thus, for r = 0, both the position and 
the maximum typically grow diffusively as ^/n for large n. 
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FIG. 1. (a) The allowed lattice points for the walker to move in the (x,m) plane. The walker is restricted to move only along the horizontal 
(constant m) lines except when it is on the x = m line. When on the x = m line, the walker can either move to a level up on the same line 
(m.m) —!> (ra+ l,m+ 1) with probability 1/2 or move to the left site while staying on the same level {m,m) —>■ (m — l,m) with probability 
1 /2. (b) The comb lattice structure after the transformation y = m~x. In the (y,m) plane, the walker is restricted to move along the vertical 
lines except when y = 0. From a point (0, m) the walker can go either to (0, ra + 1) with probability 1 /2 or to (1, m) with probability 1/2. 


When the resetting to the maximum is switched on (r > 0), 
the walker feels a dynamical bias towards the maximum. 
Hence, one expects that both x{n) and m{n) will grow faster 
than pure diffusion for large n. The question is how much 
faster? We will see that for r > 0 (strictly), both m{n) and 
x(n) grow linearly with n for large n with the same speed. In 
addition, the variance of both m{n) andx(n) grow diffusively 
for large n with the same diffusion coefficient. This suggests 
that for all r > 0, the position latches on to the maximum 
and indeed, we show that the difference variable m{n) —x{n) 
approaches a stationary distribution as n —> oo for all r > 0. 
Our exact results are summarized below. 

Statistics of the maximum m{n): The average maximum, for 
large n, behaves as 


where D{r), for r > 0, is given by 


D{r) = 


{l-r)r 


^r(2 — r) r — 2r^ + y/r(2 — r) 

(2-2r- 5r^ + 3r^) + (2-r-+ 2r^)y/7(2 


( 10 ) 


r) 


The behavior of v(r) and D(r), vs. r, are shown in Fig. 
Note that as r —>■ 0, D(r) —>■1/2 7 ^ D,„(0) = (1 —21%). Thus, 
there is a discontinuity in the variance as r —>^ 0. These results 
clearly indicate that r = 0 is a singular/critical point. Indeed, 
we find that near the critical point r — 0, there is a scaling 
regime. Taking r —?> 0, n —> oo, but keeping the product rn 
fixed, we find that the mean and the variance exhibit the fol¬ 
lowing scaling behavior 


{m{n)) ~ Y — for r = 0, 

~v{r)n for r > 0, 

where the speed v(r) is given by 

r(l-r) 


<r) = 


r — 2r^ + y/ r{2 — r) 


(4) 

(5) 

( 6 ) 


The speed vanishes as v(r) « y/ r/2 as r —> 0. The variance of 
m{n) grows diffusively for large n. 


= ('«^(«)) - ('«(«))^ - Dm{r)n, 


(7) 


where the diffusion coefficient (the subscript m in Dm{r) de¬ 
notes the random variable m) is given by 


0-^) 

for r = 0, 

(8) 

Dir) 

for r > 0, 

(9) 


{m{n)) ^ y/nfm{rn), (11) 

a^^nFm{rn), ( 12 ) 

where the two scaling functions fm{y) and F^iy) have non¬ 
trivial expressions 


fm{y) — 

Fm{y) ~ ^ F 2 



fmiy)^ 


erf(v6') + 



(13) 

(14) 


where erf(z) = -^ Jo ^ du, is the error function. The scal¬ 
ing function fmij) have the asymptotic behaviors: fmiy) 
y/2l% F y/tl%{y 13) asy^O, and/„,(y) -- 
as y —>■ oo. Consequently, Fmiy) —>■ (1 — 2-1%) as y —0, and 
FmJy') —>■ 1/2 as y —>■ oo. These two limiting behaviors of 
the scaling functions (for the mean and the variance) then 
smoothly interpolate between r = 0 (strictly) and r > 0 (and 
n oo). For any small but fixed r, there is a crossover time 
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FIG. 2. (Color online) v(r) and D(r) as a function of r (in blue and 
red respectively). The dashed lines show their limiting values, 1 /3 
and 8/27 respectively, for r —>■ 1. 


The scaling functions have the exact expressions 


fx{y) 

Fx{y) 


1 

T 
2 


(y-0 erf(v^) + 


\—e 


-fKy)- 



( 22 ) 

(23) 


The scaling function fx{y) has the asymptotic behavior: 
fx{y) ~ - (2/15)v/^/ as y 0 and 

/r(y) ~ \/y/2 — 1 /v/8y as y —oo. As a result, F^ij) —>■ 1 as 
y —>■ 0 and Fx{y) ^ 1/2 as y —>■ oo. These scaling functions 
then smoothly interpolate between the critical (r = 0) and 
off-critical (r > 0) growth of the mean and the variance. 


Statistics of the difference y(n) = m(n) —x[n): We show that 
the probability Qy{y,n) that the position at the n-th step is at 
a distance y away from the global maximum, has the large 
deviation form 


gy(y = wn,n)~exp[-«//(w)], (24) 


n*(r) 1/r, such that for n < n*{r), the mean and the vari¬ 

ance grow with n simply as a random walk (r = 0): {m{n)) ^ 
y/lnpK and (o;^,) ~ (1 — 2/;r)n. However, for n > n*(r), the 
walker starts sensing the presence of a finite resetting rate r 
and crosses over to a new behavior where the mean max¬ 
imum grows linearly with n, {m{n)) ~ v{r)n and the vari¬ 
ance grows diffusively, ^ D{r)n with the diffusion con¬ 
stant D{r) given in Eq. (10 1 . 


Statistics of the position x{ny. We find that the mean position 
behaves as 


{x{n))=0 for r = 0, (15) 

~v{r)n for r > 0, (16) 

where the speed v(r), for r > 0, is the same as that of the max¬ 

imum given in Eq. (j^. The variance of x{n) grows diffusively 
for large n, 

Gx = {x^{n))-{x{n))^ :^Dx{r)n, (17) 

with the diffusion coefficient (the subscript x denotes the ran¬ 
dom variable x) 

Dx{r) = l for r = 0, (18) 

= D{r) for r > 0, (19) 

where the diffusion coefficient D{r) is the same as that of the 
maximum and is given in Eq. ([T§. As in the case of the max¬ 
imum, the diffusion coefficient D(r —0) = 1/2 Djc(O) = 1 
undergoes a discontinuous jump at the critical point r = 0. 
Similar to m{n), the behavior of the mean and the variance of 
x{n) exhibit a scaling behavior in the scaling regime r ^ 0, 
n —y oo with the product rn fixed 

{x{n)) x/nfx{rn), (20) 

Gx ^ nFx{rn). (21) 


where the large deviation function is given by 

' 14 - 7 ^ 2 ^ 


win 


H{w) = < 


1-r 


w 1 -F w v/l — 

In - 


f 2 1 — w 


1-r 


for w < w*, 


for w > w*, 


(25) 


with w* = ^ r(2 — r). This result shows that for a given n, the 
probability Qy{y,n) becomes independent of n fory < w*n: 


2y(y,n) ~exp 


-In 




1-r 


(26) 


while for y > w*n, the distribution Qy{y,n) is still time- 
dependent. In other words, the distribution of y becomes 
stationary on a larger and larger length scale y*(n) = w*n 
that grows linearly with time n. Moreover, the rate function 
H{w) is weakly non-analytic at w = w*: the second derivative 
H"{w) is discontinuous at w = w*. This signals a dynamical 
phase transition, similar to the one observed in the temporal 
evolution of the distribution of position of a Brownian motion 
with resetting to its initial position im. 


III. THE DERIVATION USING GENERATING 
FUNCTIONS 

In this section we outline the derivation of our results. We 
start with the dynamics of the two basic observables x{n) 
and m{n) given in Eqs. (j^ and ([^. Since, we have m > 0 
and X < m, it is convenient to define the difference variable 
y = m—x, where y > 0. In terms of y, the dynamics in Eqs. 0 
and Q get translated into the equivalent forms [see Eig.[T](b)]: 
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ify > 0 


{y — I,m) with probability (1 — r)/2, 
{y, m) —(y + 1, m) with probability (1 — r) /2, 
(0, m) with probability r, 


(27) 


and if y = 0 
(y,m) - 


(0, m + 1) with probability 1 /2, 
(1, m) with probability 1 /2. 


(28) 


Let P{x,m,n) and Q{y,m,n) denote the joint probability 
distribution of (x, m) and (y, m) respectively, at the «-th time 
step. Evidently, P{x, m, n) = Q{m—x, m, n). Using the dynam¬ 
ics in Eqs. ( |27] l and ( |28| l, it is easy to write down the master 
equation for Q{y,m,n) as 


Q{y,m,n) = 


for y > 0, and 


1-r 




Q{y- 


+ ^^^Q{y+l,m,n-l), (29) 


Q{0,m,n) = 1) -f ^0(0,m- l,n- 1) 

+rf^Qiy,m,n-l), (30) 
.v=l 

with the initial condition Q{y,m,Q) — 5yfl5,nfl, and the bound¬ 
ary conditions Q{y = 0 and Q{y,m —>■ oo,n) = 0. 

To solve the set of linear equations ( [29] l and (30 1 , it is natu¬ 
ral to define the generating function 


G{s,z,X) = Q{y,in,n)s^z"'X’'. 

y=0m=0fi=0 


(31) 


Evidently, G{\,z,X) is the generating function for the prob¬ 
ability distribution of the global maximum position and 
G(s, 1,A) is the generating function for the probability dis¬ 
tribution of the position. Moreover, G(l, 1,A) must be equal 
to (1 — A)^* as demanded by the normalization of the proba¬ 
bility. After straightforward algebra, it follows that G(s,z,A) 
satisfies 


G(i,z,A) 


a / 1 

‘-2 


1 -f ( rs---hz —2rj —E(z,A) -|-rAG(l,z,A), (32) 


where 


a = (1 — r)A 


(33) 


F{z^)=Y^Y^Q{0,m,n)z”^X\ (34) 

«=0m=0 


The expression for G(l,z,A) can be obtained by setting s = 1 
in Eq. ([32]i as 


1 


G(l,z,A) = ^ 


l-(l-z)-E(z,A) 


(35) 


The normalization condition G(1,1,A) = (1 — A)^' is im¬ 
mediately checked from above. Substituting G(l,z,A) in 
Eq. ( |3^ , after some algebra, we get 

G(s,z,A)^(s+-s)(s-s_) = 


l-a 

1-A 


1-[zi(i)-z]^E(z,A)^ , (36) 


where 


and 


s± = 


1 ± v/l — 


Zl(i) = 


1-A 
I — a 


1-r ^ rA 

S I — A 


(37) 


(38) 


In Eq. (36 1 , the function F{z,X) is still undetermined and 


has to be determined self-consistently. To proceed, we note 
from Eq. (36 1 that G(i,z,A) has two poles at i = respec¬ 
tively. Therefore, inverting with respect to s, gives the form 


oo oo A R 

E Ee(A'«,«)z"’A« = - + - 

A=0m=0 ■*+ ■*- 


(39) 


where A and B are the residues at the two poles. However, 


from Eq. (37 1 , we notice that > 1 and < 1. Hence, 
the (1 /s-)' term in the above expression diverges when y — 
oo, which is inconsistent with the boundary condition Q{y 
oo,m,n) = 0. The only way to prevent this blow up is that 
B must necessarily vanish, which implies that the right hand 
side of Eq. ( |3^ vanishes for s = S-. We note that this method 
of determining the self-consistency condition via the ‘pole¬ 
cancelling’ mechanism was used before in other contexts ll27l - 
|29l- This condition then determines F{z,X) as 


A , ,, 1 

-E(z,A) =-, 

2 zo-z 


where zo = zi(i-)- (40) 


Substituting F{z,X) in Eq. (351, gives the generating function 
of the maximum as 


E E^m(m,n)z'”A"=G(l,z,A) = - 


1 zo-l 


m=0n=0 


-X zo-z 


(41) 


The normalization condition G(1,1,A) = (1—A)^* is imme¬ 
diately checked. It is easy to invert G(l,z, A) with respect to z 
exactly, which gives 


E A’m(»J,n)A" = 

n=Q 


ZO-l 1 


1-A 7™+' ’ 


(42) 


where Pyn{m,n) is the probability of having the maximum at 
m inn steps. 


and 
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On the other hand, substituting F(z,A) in Eq. (36i yields 
the full generating function as 


G(s,z,A) = 


2s+ 1 — r(l — ss_) 


(43) 


« (i+-i)(zo-z)’ 
where we have used = 1 and 
zo-ziC*) = (s-s-)s^'s+ [l -r(l -ss_)] (1 - A)(l 
Using zo = Zi (■?-), it can be shown that 

(zo- 1)(1 -S-) -A)[l -K1 -*-)]■ (44) 

Therefore, the normalization condition G(1,1,A) = (1—A)^* 
is checked from the above expression. Finally, it is useful to 
rewrite the above expression as 


G(i,z,A) = 


■^+ 

s+ — s 


1 

zo-z’ 


where we recall that Eqs. ([^ and 


zo = 


(1-A) 

(1-a) 


(I-?-) . rX 

- - -+2r + 


(1-A) 


(45) 


(46) 


using Eq. (38 i in zo = zi(s-). The explicit expression of the 
generating function G(s,z,A) in Eq. (45 i is the central result 
of this paper and is valid for arbitrary 0 < r < 1. 

Let us first check a few immediate consequences of this 
result in Eq. One can easily invert Eq. ( [45] l with respect 
to s and z to give 


Y, Q(.y^m,n)X" = - i/ - r5v.i 

n 


«=0 


-(m+1) 


Now, using P{x,m,n) = Q{m — x,m,n) we get 

Y P{x,m,n)X'^ = — [4 (5+zo)-'" - r5,,,„Zo'”] ■ 
„to azo 


(47) 


(48) 


From Eq. ( |48| ), one can derive the marginal distribution of 
the position by summing over m. Note that for positive x, the 
sum over m goes from x to In contrast, for negative x, the 
sum goes from zero to This yields 


Yp.MX'' = - 
“o a(s+zo-l) 


[(1-r)i+zo + r]zo^^^*^ forx>0, 


J^+i) 


(49) 


forx < 0. 


For completeness, we also compute the generating function 
with respect to both x and n, which gives 


a (zo-s)(s-s-) 


(50) 


As mentioned above, the expression Eq. ( |45] ) for the full 
generating function G(s,z,A), valid for arbitrary reset proba¬ 
bility 0 < r < 1, is the main central result of our paper. Several 
asymptotic results for the statistics of the two random vari¬ 
ables m and x can then be derived by analyzing this expres¬ 
sion in Eq. ( |45| ) in different limits, which we present in the 
next few sections. 


IV. THE r = 0 CASE 

Let us first check the case without resetting, i.e., we set 
r = 0. This is the standard one dimensional random walk and 
the results for the maximum and the position of the walker are 
well known. However, we reproduce it here as a check as well 
as for the sake of completeness. 

For r = 0, a = A, and 


1 -(- 4l — A^ 
Zo = ■S+ = --s-■ 


(52) 


Putting X = e P, and taking p —> 0 limit, we have zo ~ 
1 -1- y/Jp, (zo - ly (1 - A) ^ \/ (2/p) and Zq -> exp(my2p) 
as m —oo keeping niy/p fixed. Therefore, in this limit, 
from Eq. ( |4^ 


YPm{m,n)e 



(53) 


Inverting the Laplace transform gives. 


A’m ('«,«) 



where m>0. 


(54) 


From this distribution, it is easy to compute all the moments 
for large n. For instance, the mean and the variance of the 
maximum grow asymptotically as 


»*(«)) - \ — and (7^ ~ 1 - - n 


(55) 


i stated respectively in Eqs. (j^ and ([^. 
Similarly, from Eq. (|49]l, we get. 


YPx{x,n)e 


-pn 


V2p 


exp 


(-|x|4^), 


(56) 


which, after Laplace inversion, gives the expected Gaussian 
distribution (since a random walk, for large n, converges to a 
Brownian motion) 


Px{x,n) 



(57) 


Setting s = 1 above and using Eq. ( |44) i, it is easy to verify the 
normalization condition 

Y f;Px(u«)A” = (l-A)-'. (51) 

A—— oo ^=0 


Thus the mean and the variance of the position behave as 
(x(n)) = 0 and (7^ ~ , 

as stated respectively in Eqs. ([^ and ([g. 


( 58 ) 
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For the joint distribution of the maximum and the position, 
we get from Eq. (|48l) 






' 2 exp 


«=o 


— {2m — x)\/Yp 


(59) 


(2m —x)'^ 
2n 


(60) 


which, after Laplace inversion, gives 
{2m —x) [2 

where x <m and m>Q. 

Let us also verify Eq. (50 1 using the exact result of the ran¬ 
dom walk 

Px{x,n) = , when n-fx is even, (61) 


and zero otherwise. Therefore, 


^ P^(x,n)s^ = Y^P{2m-n,n] 


„2m—n 


m=0 


1 / 1 


2 1 i 


Consequently, 


^ Px(^,«K = 

n=0 x=-n 


, >■ f 1 


2s 


1 


X {s+ — s)(s — S-) ’ 


(62) 


which is same as Eq. (50 1 for r = 0. 

Thus, for the case r = 0, when there is no resetting to the 
maximum, we have recovered the results for the usual random 
walk. 


V. THE r 1 LIMIT 

In the other extreme limit of r —>^ 1, we have zo = (2- 
X'^)/X. Therefore, from Eq. (42 1 


Y^Pm{m,n)X" = 

n=0 


2 + X 

2^ 


2-A2 


(63) 


From the series expansion of the above expression with re¬ 
spect to A, it is easily verified that P^{m,0) = 5„t.o- Moreover, 


PmiQ,n)=2 


_ («+!) 


Wzl)!+V2i±H)! 


(64) 


is the probability that the walker has not crossed the origin up 
to the step n. This is the case where the walker goes to the site 
X = — 1 (with probability 1/2) and comes back to the origin 
(with probability one) at alternate time steps. In general (for 
r ^ 1), 




A 

2-A2 


(65) 


For large m and n, with m = wn, the integral can be evaluated 
using saddle-point approximation, which to the leading order 
gives 


Pm(m = wn,n) ^ exp[—n5(w)] 
where 5(w) = 5(tv,A*) with 

A 


5(w,A) = InA — win 


2-A2 


( 66 ) 


(67) 


The saddle point A* is obtained by solving the equation 
dxS{w,X)\x* = 0, as [A*(w)]2 = 2(1 — w)/(l -Fw). Substitut¬ 
ing this in the above equation we obtain the large deviation 
function as 




4w 

-win-—-. (68) 
(1 +w) 


This large deviation function has a maximum at w = 1 /3, and 
near this, one gets S{w) = (27/16)(w— 1/3)2, which implies 
the Gaussian form 


~exp - 


In fact, Eq. 


(m — nj3y 


2Dn 


withD=^. (69) 


can be inverted exactly, which gives 

if (n — m) is even, 


(n+ra) 

2 2 - 


Pm{m,n) = 


n+m 

2 

m 


n+m+l) 

m 


(70) 


2 ^ 2 ^ ^ \ if (fi — m) is odd. 

For large n and m, using the Stirling’s approximation in 


Eq. (70 1 , one can recover the large deviation function given 
by Eq. (68]l. 

Let us now look at the probability distribution of the posi¬ 
tion using Eq. (|4gi. Eor r —> 1 we have a —> 0 and 00 
with 2. Therefore, from Eq. ( [49| , it is clear that for 

negative x, we get nonzero probability only for x = — 1, which 
reads 


f Px(-l,«)A" = - = ;^. 
to ^0 2-A2 

Inverting this with respect to A gives 

7^x(-l,«) =- 


(71) 


(72) 


This result can be understood, as this is the case where the 
walker goes to the site x = — 1 (with probability 1 /2) at odd 
time steps and comes back to the origin (with probability one) 
at even time steps. The also implies that P{—l,n) is nonzero 
only when the maximum remains zero. Indeed from Eq. ( [48| , 
for r ^ 1 we get l)„F(-l,0,w)A" = Zq ' = L«A’x(-l,«)A". 


For X > 0, Eq. (|49|l we get 




A 

2-A2 


(73) 
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Although the exact form differs from Eq. ( [6?| , it is evident that 
Px{x,n) and Pm{m,n) have the same large deviation function, 
i.e., 


Px{x = wn,n) ~ exp[—n5(w)], 


(74) 


with S{w) given by Eq. ( 681 . 


VI. THE CASE 0 < r < 1 

Eor general resetting probability 0 < r < 1, it is a bit cum¬ 
bersome to find the exact large deviation functions associated 
with the probabilities P^{m,n) and Px{x,n). However, one ex¬ 
pects (as in the cases of r = 0,1) the typical fluctuations near 
the mean to be governed by Gaussian distributions. The gen¬ 
erating functions for the mean and the variance for the max¬ 
imum and the current position are obtained from the gener¬ 
ating functions of their distributions, given by Eq. •ED and 
Eq. ( [soi l respectively, simply by taking derivatives. Erom their 
respective generating functions, we then derive, for arbitrary 
but fixed 0 < r < 1, the asymptotic behavior of the mean and 
variance for large n for both the maximum and the position. 
Einally, the asymptotic behavior of the the full probability dis¬ 
tribution of the difference variable y(n) = m(n) —x{n) is de¬ 
rived. These derivations are outlined in the next three subsec¬ 
tions. 


A. Statistics of the maximum m{n) 


Erom Eq. we get the generating functions of the first 
two moments of the maximum as 


n=0 


dz 


G(l,z,A) 


Z=1 


1 1 
1 - A zo - 1 ’ 


(75) 


and 


H=0 


dz^ dz 


G(l,z,A) 


Z=\ 


1-A 


(zo-1)^ zo-l_ 


(76) 


respectively. Using the expression of zo = from 

Eq. ([3^, it is easy to check that 


1 C{r,X) 

zo-1 1-A ’ 


(77) 


where 


C(r,A) = 


j(l — a) 


(1 — a)(l — 2r) -fv/l — 


(78) 


with a = (1 — r)A. Inverting Eqs. (75 i and (76 1 with respect 
to A using Cauchy’s formula, we get 


(m(n)) 


1 f dX C(r,A) 
iTZi Jo A"+l (1 - A)2 


and 


\n)) = 


1 , 

[ 

r2C2(r,A) C(r,A) 1 

27ti 7 

0 A"+i 

[(1_A)3 (1-A)2j 


(80) 


respectively, where the integral is along a counterclockwise 
closed contour around the origin in the complex A plane. 

A hint that r = 0 is a special case and is different from the 
r > 0 case can be already seen at this level. Eor r = 0, a = A 
and from Eq. (78 1 one gets C(0, A) = (1 — A)A [(1 — A) -f 

\/l — A^] '. Therefore, for the r = 0 case, the powers of (1 — 
A) in the denominators of the above expressions, reduce by 
one, and hence this case must be treated separately from the 
r > 0 case. In the following we will only consider the case 
r > 0. 


Evaluating the contour integrals in Eqs. (79 1 and (80i ex¬ 
plicitly for all n looks cumbersome. However, the asymp¬ 
totic behavior for large n can be derived by separating out the 
contributions to the contour integrals arising from the pole 
at A = 1. Note that C(r,A) also has two branch points at 
A = ±(1 — r)^'. Therefore, fo — + [contributions from 

the integrals around the branch cuts from —°o to —(1 — 
and from (1 —r)^* tooo]. Using the residue theorem and com¬ 
puting the residue at A = 1 we get 


1 f dX C(r,X) „ , / X 

-(- [branch cuts contributions], 


(81) 


and 


1 / dX 2C^(r,X) , ,,, X 

^JoA^+i (1-A)3 -(«+l)(« + 2)Co(^) 

+ 4(n-f l)Co(r)Ci(r) +2C?(r) +2Co(r)C2(r) 

-b [branch cuts contributions], (82) 


where 




dX^ 


C(r,A) 


A=1 


d'' 

dX^ 


C(r, 1-A) 


A=0 


■ (83) 


These coefficients can be calculated explicitly using Eq. Gi- 
Eor example, the first two coefficients are given explicitly as 


Co(r) 


r{l-r) 

r (1 — 2r) + \/2r — 


(84) 


and 


Ci(r) = 


r(l—r) 1 — 3r-|- -|-r(2r — \ )\/2r — 


^/2r — r^ (^r — 2r^ -\- \/2r— 


(85) 


One can also show that the branch cuts contributions in Eqs. 
and ( |82| ) go to zero exponentially fast as n —> oo. 

Using the results from Eqs. ( |8 T] i and ( |8^ in Eqs. ( |79| l and 
([80]), we then obtain, for large n, the mean 


{m{n)) ~ v{r)n, with v(r) = Co(r) 


r{l — 2r) -b 'J2r — 

( 86 ) 


(79) 
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The speed v(r), as a function of r, is plotted in Fig. Simi¬ 
larly, the variance grows linearly for large n 


= ('«^) - ('«)^ - D,„{r)n, (87) 


with 


D„,(r) =D(r) = Co"(r) +2Co(r)Ci(r)+Co(r). ( 88 ) 


Using Co(r) and Ci(r) from Eqs. (84i and (851 gives the ex¬ 
plicit expression for the diffusion coefficient D(r), for r > 0, 
as given by Eq. (10 1 . A plot of D(r) vs. r is provided in Fig.|^ 


B. Statistics of the position x(n) 

Similarly, to compute the mean and variance of the position 
x(n), we take derivatives with respect to s in Eq. (50l at s = 1 
to get 




n=0 

and 


1—r(l—i_) 

(l-s_)(zo-l)^ (l-s_)^(zo-l). 


(89) 




n=0 


2 [l-r(l-s)] 


(l-._)(zo-ir 


S- (1 +S-) 


1—r(l—s_) — 3s_ 

where we recall Eqs. ( |3T| , ( |T7| i and ( |46| ). 

It is useful to define the following quantities 


(90) 


Ai(r,A) = 
AiinA) = 
A3(r,A) = 
A4(r,A) = 


2 [l-r(l-^_)] 
a(l—s_) 

2s- 


j(l — ’ 


2 [l-r(l-5_) -35-] 
a(l—5-)^ 

2s_ (1 -|-s_) 

a (1 — i_) 


3 • 


(91) 

(92) 

(93) 

(94) 


We note that Ai(r,A)C(r,A) = 1 and [ 2 A 2 (r,A) -f 
A 3 (r, A)]C(r, A) = 1. Using these definitions and Eq. ( [77| , we 
invert the above generating function with respect to A and get 




C(r,A) A 2 (r,A)C(r,A) 


(1-A)2 


( 1 -A) 


, (95) 


2C2(r,A) A 3 (r, A )c2(r, A) 

(i-A)3^ [i-xy 

A 4 (r, A )C(r, A) 


-f ■ 


( 1 -A) 


(96) 


Einally, evaluating these integrals using residue theorem, we 
get 


and 

(x^(n)) = (n-f l)(n-|-2)Co 

-I- (n-f 1) [4Co(r)Ci(r) +A^fi{r)Cl{r)] 
+ 2 Cf(r)+ 2 Co(r)C 2 (r)+A 3 .i(r)C^(r) 
-I- 2 A 3 ,o(r)Co(r)Ci (r) -f A 4 .o(r)Co(r) -f . 
where C„ (r) is defined by Eq. (| 8 ^ and 


(98) 




A/.„(r) = (-l)"^A,(a) 


X = \ 


d’‘ 

= ^A;(r,l-A) 

<9A" ^ ' 


x=o 

(99) 

Einally, collecting the leading term for large n, we obtain 
the following asymptotic results for the mean and the vari¬ 
ance. Eor large n, and for r > 0, the mean 


(x(n)) ~ v{r)n, 


( 100 ) 


where the speed v(r) = Co(r) turns out to be exactly the same 
as in case of maximum, with explicit expression given in Eq. 
( | 86 | ). Similarly, the variance behaves, asymptotically for large 
n and for all r > 0 , as 

0-2 = (x^(n)) - {x{n))^ Dj,{r)n, 

with D^(r)=D(r)=C^(r)+2Co(r)Ci(r)+Co(r). (101) 

Hence, the first and the second moments of the two variables 
m{n) and x{n) grow identically for large n. It is then natu¬ 
ral to ask how the difference variable y{n) = m{n) —x(n) is 
distributed and we address this in the next subsection 


C. Asymptotic distribution of the difference y{n) = m{n 


K{n) 


Let Qy(y,n) be the probability that the location at the n-th 
step is at a distance y away from the global maximum. Clearly, 
Qy{y,n) — TZ=oQiyAn,n). Therefore, fromEq. (471 we get 




n=Q 


1 

-I zo-L 


( 102 ) 


Let us first look at the value y = 0, for which Eq. ( fT02| ) 
reads 


^(2y(0,n)A" = 


1 


n=0 


_ 2C(r,A) 
A ZQ — 1 A(1 — A) 


(103) 


where we have used Eq. ( 0 . To extract the large n behav¬ 
ior of Qy {0,n), we need to investigate the right hand side of 
Eq. ( |103| ) in the limit A —?► 1. As A —>■ 1, C(r,A) —>■ C(r, 1) = 
Co(r) given in Eq. (84i. Hence, the right hand side of 
Eq. ( |103 1 behaves as 2Co(r)/(l — A) as A —>■ 1. This clearly 
indicates that Qy{Q,n) becomes independent of n in the limit 
n —y 00, given by 


Qy{0,n -)■ 00 ) -)■ 2Co{r). 


(104) 


Eor y > 0, inverting the generating function Eq. (|102 1 


we 


get 


, . 1 f dX 2 


(x(n)) = (n + 1 )Co (r)-|-Cl (r)-A 2 ,o(r)Co(r)(97) 


(105) 
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Now, using explicit expressions, and changing the integral Finally, combining these two regimes we obtain the large de- 
over X to that over a = (1 — r)X, we get viation behavior 


2{l-a) 

(1 — a)(2ao— 1) + Vl — 


/■ao 

a 

V (3 / 

. 1 -f Vl — a^. 


where ao = I — r. The integrand has a simple pole at a = ao 
and branch points at a = ±1. The contour of integration 
around zero can be split into two vertical contours: one that 
goes from +i°o to —i°° through the left of the origin and an¬ 
other that goes from —i°o to +i°° through the right of the ori¬ 
gin. The left contour subsequently can be wrapped around 
the branch cut from a = — 1 to —oo. The contribution from 
this contour is subdominant and the main contribution comes 
from the contour on the right for large n. It is useful to express 
the integral as 


, , 1 / da 


Qy{y = wn,n) 


1 

2m 


|■{0+)+ioo 
I (0-p)—/CO 


da 


gja) 

(ao-a) 


exp[—n//(w,a)], 


(107) 

so that for large n we can use the saddle point approximation 
method. Here 


H{w,a) = ln(a/flo) — win 


a 

1 -f Vl — 


(108) 


and 


8{a) 


2(1-a) 

(1 — a) (2ao — 1) + Vl — 


(109) 


The saddle point a* is obtained by solving the condition 
daH{w,a)\a* = 0, which gives a* = a/ 1 — w^. Note that the 
integrand of Eq. (|107|i h as a simple pole at a = ao = 1 — r. 
For w <w* = \Jr(2 — r), we have a* > ao, and therefore, the 
contribution to the above integral comes from both the pole 
and the saddle point. However, the contribution from the pole 
is larger than that from the saddle point. Therefore, for large 
Qy(y = wn,n) ~ exp[—n//i (w)], where the large deviation 
function is given by 


Qy{y = wn,n) ~ exp[-nH{w )], 
with the rate function given by 


H\{w) = win 


H{w) = < 


\ + ^fH2^ 


1-r 


( 112 ) 

for w < w*, 
for w > w*, 

^_ (113) 

with w* = ^r{2 — r). 

As discussed in Section 2, this result indicates that for a 
given large n, Qy{y, n) becomes independent of n for y < w*n 
and is still n-dependent for y > w* n, signalling a dynamical 
phase transition. The rate function H[w), plotted in Fig. is 
weakly singular at the critical point w = w* where both H[w) 
and H'{w) are continuous, but the second derivative H"{w) 


, . W 1-fW Vl-w2 

H 2 {w) = -In---fin—- 

2 1 — w 1 — r 


is discontinuous: w* ) = 0, while //"(w 


w 


^) = 


1 /(1 — r)^. Such a second order dynamical phase transition 
was also observed recently in the time evolution of the dis¬ 
tribution of position of a Brownian motion in one dimension 
with resetting to its initial position CD. 

We finish this subsection by making a couple of interesting 
observations. In the limit r —>^ 0 we have w* = 0. Therefore, 
the large deviation form is given by 7/2 (w) with r = 0. Ex¬ 
panding in Taylor series, for small w we get H 2 {w) ~ w^/2, 
which gives the Gaussian form Qy{y,n) ^ exp[—y^/(2n)]. 
This is consistent with the result obtained by substituting 


X = m — y in Eq. (60 1 and integrating over m from 0 to 00 . 


On the other hand, in the limit r ^ 1, we have w* = 1. There¬ 
fore, the large deviation form is given by Hi (w) with r — 1. 
We note that for r = 1, the large deviation function is — 00 , ex¬ 
cept for the case w = 0. This is because, y can take only two 
values, namely, 0 and 1, for the case r = 1. Using the result of 
Eq. (104i, we get Qy{0,n —> 00 ) —2Co(l) =2/3. For y = 1, 
Eq. (102 1 gives 


n=0 


2 

{as+)izo -1) ■ 


(114) 


Hi (w) = H(w,ao) = win 


1 + 7K2^ 

1-r^" 


Now, in the limit r —> 1, we get as+ —)• 2 and (zo — 1) —> 
(110) A^*(l — A)(2-f A). Using these, and inverting the above 

equation we get 


On the other hand, for w > w*, we have a* < oq. Therefore, 
evaluating the above integral using saddle-point approxima¬ 
tion gives Qy{y = wn,n) ~ exp[—n// 2 (w)] where the large de¬ 
viation function is given by 


H 2 {w) =//(w,a*) = ^In I + ln 
2 I —w 




( 111 ) 


2y(l,«) 


1 

3 



as n —>• 00 . 


(115) 


Finally, following the method used in Ref. OOl in a differ¬ 
ent context, we can also write down a more complete asymp¬ 
totic form of Q{y,n) for large n as 


1-r 
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Qy{y = wn,n) : 


g-nH2(w) 

2y/nn 


K[w) — 


sgn(w* - w)g-i 




sgn(w* — w) 


erfc(^-\/n[H2{w) — Hi (w)]) — 0(w* — w) 

(116) 



FIG. 3. (Color online) Q(y,n) against the scaled variable w = yjn 
for n = 100 and r = 0.01. The (blue) points are obtained from nu¬ 
merical simulation. The solid (red) line represents the asymptotic 
form given by Eq. l |l 16| ). The dashed lines (magenta and black) plot 

with the 


the large deviation forms Qy{y = wn,n) ~ exp[— 
large deviation functions Hi (w) and H 2 {w) given by Eq. 1110 1 and 


Eq.fTTT) respectively. The vertical dotted line marks the position of 
: 0.141. 


where g_i = —g(ao) andA'(w) = V2wg{a*)/{ao — a*). This 
pre-asymptotic form is particularly useful to compare to the 
results of simulation. Indeed, Figure compares this form 
with the numerical simulation results. The agreement is ex¬ 
cellent. 


VII. THE LIMIT r ^ 0 AND THE ASSOCIATED SCALING 
FUNCTIONS 


From the results presented in Section IV for r = 0 and in 
Section VI for r > 0, we see that the limit r—>-0 (after taking 
the large n limit) is not the same as r = 0, for the statistics of 
both the maximum and the position. In other words, the two 
limits limr^o lim„^.,o do not commute. This indicates that 
r = 0 is a singular or a ‘critical’ point. For finite but large n, 
there should then be a smooth crossover function interpolat¬ 
ing between these two limits. In this section, these crossover 
scalings functions are derived analytically and compared to 
numerical simulations. 


A. Scaling functions associated with the maximum m(n) 


We first consider the crossover scaling functions (near r —>^ 
0) associated with the mean and the variance of the maximum 
m(n). Let us first focus on the mean. Our starting point is the 
exact generating function for the mean in Eq. ( |75] l, where we 
recall that zo = Zi (s-) is given in Eq. (381. We also remind 
the reader that 


a = (1 — r)A and = - ^1 — y/1 — . 


(117) 


Since, we want to analyze the behavior of {m{n)) for large n, 
and simultaneously r —>^ 0, we need to set X close to 1 in Eq. 
(751. We set, X = \— p where p is small. Upon inspecting 
S- in Eq. it follows that the right scaling limit is when 

p —> 0, r —0, keeping the ratio pjr fixed. In the real space, 
this limit corresponds to r —> 0, « —“o but keeping the product 
rn fixed. In this limit, the leading behavior of can be easily 
worked out to give 


s-^l-y/2{r + p) . (118) 

Substituting this leading behavior of s_ in Eq. ( |38] l, it follows 
that in the scaling limit 


zo =zi(s-) « 1 + 


VXp 

y/r + P ' 


(119) 


We substitute this leading behavior of zo on the right hand side 
(rhs) of Eq. (75 i and use A = 1 — p. In the limit p —> 0, the sum 
on the left hand side (Ihs) of Eq. ( [75] l can be approximated by 
an integral, yielding in the scaling limit 



( 120 ) 


By power counting on both sides of Eq. ( |120| i, it follows that 
(m(n)) must have the following scaling behavior 


(«(«))«V«/,„(rn), (121) 


in the appropriate scaling limit r ^ 0, n ^ oo with the product 
rn fixed. 

Substituting this scaling behavior on the Ihs of Eq. ( |120| i, 
setting pjr = s and making a change of variable rn = y, yields 
the following equation for the scaling function fmiy) 

( 122 ) 
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To invert the Laplace transform on the rhs of Eq. ( 122| i, 
reexpress 

\/l + i 


1 


+ ■ 


1 


(123) 


sVT+s 

Each term on the rhs of Eq. ( |123| l is an elementary function 
that can be easily inverted using convolution theorem. Invert¬ 
ing, we then get an exact expression for the scaling function 






^ -y 
e ’ 


(124) 


where erf(z) = ;^ /o ^ “ du- The function fm{y) has the fol¬ 
lowing asymptotic behaviors 


fmiy) ~ < 


—I- 0{y) as y —y 0, 

K 


(125) 




as y —y < 


Thus, when r = 0, using /m(0) = y/^/TT in Eq. (121 1 yields 
the asymptotic behavior of the mean, {m{n)) ~ s/lnln. In 
contrast, when r > 0, as n —y oo, the scaling argument y = m^ 
oo. Hence, using the other asymptotic behavior in Eq. ( |125| l 
as y —y oo, yields the linear growth {m{n)) ~ y/ rjln. Note 
that the speed v(r) in Eq. (|6l indeed tends to v(r) ^ 
as r —y 0. The exact scalin^unction /m(y) thus interpolates 
smoothly between these two limits. Eor any small but nonzero 
r, we thus expect that {m{n)), as a function of n, will initially 

grow as ~ \l^ (the critical behaviour at r = 0), before cross¬ 
ing over at a characteristic time n*{r) ~ 1/r to the off-critical 
linear growth, {m{n)) ~ y/r/2n. In Eigure|^ we compare the 
numerical simulation results for small values of r and show 
how they approach the analytical scaling function f,„{y) as 
r-y 0. 

Next we consider the variance of the maximum, = 
{m^[n)) — {m{n)y in the scaling limit r ^ 0, n ^ o° but keep¬ 
ing the product rn fixed. Eor this, we now need to analyze the 
second moment in Eq. ( |7^ in the scaling limit. The analysis 
proceeds more or less as in the case of the mean. We do not 
repeat this computation here and just mention the final result. 
In the scaling limit r —> 0, n —y oo with the product rn fixed, 
we find that the variance behaves as 


Ctm ~ nFm{rn), 

with the scaling function F^iy) given by 

Tm(y) = l + |-/^(y), 


(126) 


(127) 


where /m(y) is given in Eq. (124 1 . The scaling function has 


the following asymptotic behaviors 


Fmiy) 



1-as y —> 0, 

71 


(128) 



FIG. 4. (Color online) plotted vs. y = rn for different small 

values of r. As r —> 0, th e curves approach the analytical scaling 
function fmiy) in Eq. 1124 1 plotted as a solid line. 



FIG. 5. (Color online) CT^/n plotted vs. y = rn for different small 
values of r. As r —> 0, the curves approach the analytical scaling 
function Fm{y) in Eq. 1 127 1 plotted as a solid line. 


Hence, for r = 0, using F’m(O) = 1 — T in Eq. (126 1 gives 


the asymptotic behavior of the variance, — (l — 4) n, i.e., 
the result for normal diffusion without resetting. In contrast, 
for r > 0, as n —y oo, the scaling a rgum ent y —y oo. Hence, us¬ 
ing Fmiy) 1/2 as y —y oo in Eq. (126 1 gives, (7^ ~ n/2. This 


agrees perfectly with the finite r result, ~ D(r) n with Z)(r) 
given in Eq. ( [T^ , since D(r —^ 0) = 1/2. The exact scaling 
function F’m(y)Aus interpolates smoothly between these two 
limits. In Eigure]^ we compare the numerical simula tion r e- 
sults with the analytical scaling function Fm{y) in Eq. (127 1 . 


as y —> oo. 








































13 


B. Scaling functions associated with the maximum x{n) 


We now turn to the scaling behavior of the mean and the 
variance of the position x{n). We start with the mean whose 
exact generating function is given in Eq. ( [89] l. To analyze the 
scaling limit r —> 0, n —> oo while keeping the product y = rn 
fixed, we follow the same procedure as in the case of the max¬ 


imum. Setting X — p with p —^ 0, and using Eq. (1191, we 

hnd that Eq. 
integral 


reduces, in the scaling limit, to the following 


/ W«))' 

Jo 


-pn 


dn 




r + p) 


(129) 


It then indicates the following scaling behavior for the mean 
position 


(x(n)) « ^/nf^{rn), 


(130) 


where fx{y), using Eq. (129 1 , satisfies 


^ y/yfx{y)e '^dy- 


(131) 


One can again easily invert the Laplace transform in Eq. ( |131| l 
to get 


fx{y) = 


1 


y/2y 

It has the asymptotics 


y- 


rf(v6')- 


K 


(132) 


fx{y) ^ < 


+ asy^o. 




(133) 


When r — 0, using fx{0) = 0, one recovers the standard 
random walk (without resetting) result, (x(n)) = 0. In con¬ 
trast, for r > 0, when n —> oo, i.e., the product y = rn —>■ °o, 
using the large y asymptotic behavior in Eq. ( 133| l, one 
gets, (x(n)) ~ x/rjln, compatible with the linear growth with 
speed v(r) in Eq. (16 1 upon noting that v(r ^ 0) = y/TJl. The 
scaling function j^) interpolates between these two limits. 
Eigure [^demonstrates how simulation results converge to the 
analytical scaling function fx(y) in Eq. (1321 as r —> 0. 


We next consider the scaling behavior of the variance of the 
position, (7^ = (x^(n)) — (x(n))^, in the scaling limit r —> 0, 
n —> oo with the product y = rn fixed. Here we analyze the 
generating function for the second moment in Eq. ( |90| in the 
scaling limit. Since the procedure is identical as in the case of 
the maximum, we skip the details and present only the result. 
We hnd that in the scaling limit, the variance of the position 
behaves as 


(134) 



rn 


FIG. 6. (Color online) plotted vs. y = rn for different small 
values of r. As r —> 0, the curves approach the analytical scaling 
function fx{y) in Eq. l|132[l plotted as a solid line. 


with the scaling function Fx{y) given by 

-fxiy). 


T7 ( \ F , 1 - e ^ 

Px{y) = 2 + - 


y 


(135) 


where fx(y) is given in Eq. (132 1 . The scaling function Fx(y) 
has the following asymptotic behaviors 


Px(y) 



(136) 


Using these asymptotic behaviors, it is again easy to check 
that Fx(y) smoothly interpolates between the critical (r = 0) 
and the off-critical (r > 0) behavior of the variance of the po¬ 
sition, for hnite but large n. In Eigure [^ we compare the nu¬ 
merical simulation results with the analytical scaling function 
Fx(y) inEq. (|135|). 


VIII. CONCLUSION 

In conclusion, we have considered a model of random walk 
in one dimension where the walker, at each time step, resets 
to the maximum of the already visited positions with a certain 
probability r. Eor r = 0, it reduces to the standard random 
walk in one dimension. The presence of a nonzero resetting 
probability r changes drastically the asymptotic behavior of 
the walker. We hnd that on average, both the position and the 
maximum move with the same speed, and we have obtained 
an exact expression for this speed v(r). The Huctuations about 
the mean is again described by the same diffusion coefficient 
D(r) for both. We also obtain the large deviation form of the 
probabilities of hnding the walker at a distance y away from 
the maximum. The associated large deviation function shows 
a second order phase transition. 

An interesting extension would be to study the walk which 
resets either to the maximum or the minimum with an equal 


(7^ « nFx{rn) 
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rn 


FIG. 7. (Color online) O^jn plotted vs. y = rn for different small 
values of r. As r 0, the curves approach the analytical scaling 
function Fx{y) in Eq. 1 135 1 plotted as a solid line. 


probability r/2. In general, one could ask the question in 
higher dimension, where the walker resets to one of the 
boundary sites of the already visited sites. Yet another exten¬ 
sion, beyond the study of the fluctuations of the position and 
the maximum studied here, concerns the study of the search 
process with resetting to the maximum. For instance, it would 
be interesting to compute the mean first-passage time to an 
immobile or a moving target for a such a random walker sub¬ 
mitted to random resetting to the maximum. 
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